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The reconstruction of scalar-field dark energy models is studied for a general Lagrangian density 
p{(j), X), where X is a kinematic term of a scalar field (j>. We implement the coupling Q between dark 
energy and dark matter and express reconstruction equations using two observables: the Hubble 
parameter H and the matter density perturbation 5m- This allows us to determine the structure of 
corresponding theoretical Lagrangian together with the coupling Q from observations. We apply our 
formula to several forms of Lagrangian and present concrete examples of reconstruction by using 
the recent Gold dataset of supernovae measurements. This analysis includes a generalized ghost 
condensate model as a way to cross a cosmological-constant boundary even for a single-field case. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 

Observations suggest that our universe has entered a 
stage of an accelerated expansion with a redshift z ^ 1 
. This has been supported by a number of recent astro- 
physical data including supernovae (SN) la 0, Cosmic 
Microwave Background (CMB) and large-scale struc- 
ture Q. These observations show that about 70% of 
the total energy density of the present universe consists 
of dark energy responsible for an accelerated expansion 
(see Rcfs. Q for recent works). The simplest candidate of 
dark energy is a cosmological constant, but this scenario 
suffers from a severe fine-tuning problem of its energy 
scale from the viewpoint of particle physics 0|. Hence 
it is natural to pursue alternative possibilities to explain 
the origin of dark energy. 

So far a wide variety of scalar-field dark energy models 
has been proposed- includin g qu intessence , k-essence 
P, tachyons , phantoms |Tol | and ghost condensates 
III El- These scenarios are distinguished from a cos- 
mological constant because of their dynamical nature of 
the equation of state of scalar fields. A typical approach 
is to predict the evolution of the Hubble parameter theo- 
retically and to check the consistency of models by com- 
paring it with observations. An alternative approach is 
to start from observational data and to reconstruct cor- 
responding theoretical Lagrangian. The latter is more 
efficient to find out best-fit models of dark energy from 
observations. 

This reconstruction is simple for a minimally coupled 
scalar field with a potential V^cf)) In fact 

one can reconstruct the potential and the equation of 
state of the field by parametrising the Hubble parameter 
H in terms of the redshift z from the luminosity distance 
Dl{z) 01 . This method can be generalized to scalar- 
tensor theories El US El, f{R) gravityjp and also 
a dark-energy fiuid with viscosity terms |23l |. In scalar- 
tensor theories a scalar field called dilaton is coupled to 
gravity, which means that an additional function F{(j)) 
exists in front of the Ricci scalar R. If the evolution 
of matter perturbations 6m is known observationally in 
addition to the Hubble parameter H{z), one can even 



determine the function F{4)) together with the potential 
of the scalar field |19l| . 

In this paper we will provide a reconstruction program 
for a very general scalar-field Lagrangian density A"), 
where X = — (V(/))^/2 is a kinematic term. We shall con- 
sider an Einstein-Hilbert action with a scalar field <p cou- 
pled to a non-relativistic barotropic fluid (dark matter) 
with a coupling Q{(j)). This coupled quintessence sce- 
nario was proposed in Ref. p3 | as an extension of scalar- 
tensor theories with a nonminimally coupled scalar field 
m El El m El. In fact the scalar-tensor action in 
Jordan frame can be transformed to the action in Ein- 
stein frame with an explicit coupling between the scalar 
field and the non-relativistic fluid. 

The presence of the coupling (5(0) means that the 
parametrisation of the Hubble rate H{z) is not sufficient 
to determine the structure of theory. We shall use the 
equation of sub-Hubble matter perturbations Sm recently 
derived in Ref. [30| for the Lagrangian density p{(j), X) 
(see also Refs. |^ El). The coupling (5(</>) is deter- 
mined once we know H{z) and 5m{z) observationally. We 
note that this places constraints on the strength of non- 
minimal couplings when the coupling Q originates from 
scalar-tensor theories. 

The observations of Sloan Digital Sky Survey (SDSS) 
and Two degree Field (2dF) galaxy clustering 34, 3j| 
provide the information of matter perturbations. Galaxy 
clustering is proportional to matter clustering on large 
scales with a constant of proportionality called bias. The 
analysis of galaxy clustering itself does not determine the 
value of bias. Hence the evolution of matter perturba- 
tions is difficult to be known unless bias is determined 
observationally. However recent observations using the 
luminosity dependence of galaxy clustering together with 
the halo mass distribution b egan to provide good data 
for the determination of bias |33. El E3 ■ 1^ expected 
that we will be able to determine the evolution of mat- 
ter perturbations, (5„j(z), from upcoming high-precision 
observations. 

The recent SN analysis using the Gold dataset im- 
plies that the parametrisation of H{z) which crosses the 
cosmological-constant boundary {w ~ —1) shows a good 
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fit to data. This crossing to the phantom region {w < —1) 
is not possible for an ordinary minimally coupled scalar 
field [p = X — V {(/>)]. This transition can occur for 
the Lagrangian density p{(j), X) in which dp/ dX changes 
sign from positive to negative, but we require nonlinear 
terms in X to realise the w = — 1 crossing Ull ^■ 
It was shown in Ref. |39l | that such a crossing is possible 
if a (phantom) dark energy fluid is coupled to a non- 
relativistic fluid with a specific coupling. This can be 
also realised in a multi-field system with a phantom and 
an ordinary scalar In this paper we shall show a 

simple one-field model [p = —X + h{(t))X'^) crossing the 
cosmological-constant boundary and perform the recon- 
struction of such a theory. 



II. RECONSTRUCTION PROGRAM 



We start with a general action given by 



respectively. By using the present ratio flom of the mat- 



S = / d^.TV^ 



\R + p{cj,,X) 



+ S„i{(j>,gt,y), (1) 



where i? is a scalar curvature, p(0, X) is a general func- 
tion of (f) and X = — (1/2)(V(?!))^. The gravitational con- 
stant, = SttG, is set to be unity. We assume that 
the field 4> is coupled to a barotropic perfect fluid with a 
coupling Q{4>) = —1/ {pm^/~-g)SSm/S(j), where p„i is the 
energy density of the fluid. We shall use a sign notation 
(-,+,+,+)• 

In a flat Friedmann-Robertson- Walker (FRW) space- 
time with a scale factor a, the flcld equations for the 
action ^ are 

= p^m + 2Xpx-p, (2) 
2H = -p„, - Prn - 2Xpx , (3) 
Pni + 'SH{pjn + Prn) = Q{(l>)Pm4> , (4) 

where H = d/a, X ~ <p^/2, px — dp/dX and Pm is the 
pressure density of the fluid. A dot denotes a derivative 
with respect to cosmic time t. In what follows we shall 
consider the case of a non-relativistic barotropic fluid, 
i.e., Pm = 0. Eq. 10} is written in an integrated form 



where 



I{<p) = exp 



(0) 



'00 / 



(5) 



(6) 



Here Pm\ ap and (f)o arc the present values of the en- 
ergy density pm, the scale factor a and the scalar field 0, 



^ I thank Alexander Vikman for pointing out that the w = —1 
crossing is hard to be realised only in the presence of the terms 
linear in X 1381. 



ter fluid and the Hubble parameter Hq, p[^ is given by 
3HqQo„i- We also define the redshift parameter z, 
z = ao/a. Then the energy density pm is written 



M 

pm 

as 1 



as 



(7) 



One has /(0) = 1 in the absence of the coupling Q{4)). 
In this case one can reconstruct the structure of theory 
by using Eqs. (0) and (0) if the evolution of the 
Hubble parameter is known from observations. This was 
actually carried out for an ordinary scalar field with a 
Lagrangian density: p = X-V{(t)) [ll[lE[l3|. When the 
field (j) is coupled to dark matter, we require additional 
information to determine the strength of the coupling. 
We shall use the equation of matter density perturbations 
for this purpose as in the case of scalar-tensor theories 

MM. 



sub- 



The evolution equation for the matter den sity contrast, 
5m = 5pm/ Prrn was recently derived in Ref. |33- On i 
Hubble scales this is given by 



2Q2(^) 



2H + Q{cp)^ 5m -T^ 1 + -^-^ PmSm^O. 

J 2 [ Px \ 

(8) 

This is a very general equation which holds for any scalar- 
field Lagrangian density ^(0, X) and also for the case in 
which the coupling Q depends on the field 4>. Eq. © can 
be solved analytically for scaling solutions [s^- In fact 
it was shown that matter perturbations are suppressed 
for a phantom {px < 0) whereas they are not for an 
ordinary field {px > 0). Generally we need to specify 
the Lagrangian density p{X, (j>) in order to solve Eq. ((SJ. 
We note that the gravitational potential $ is related with 
the matter perturbation 5m through the relation $ ~ 
— {3a'^H^/2k^)^m5m on sub-Hubblc scales [s^ [k is a 
comoving wavenumber). 

One can rewrite the equations Q , and JHl by using 
a dimensionless quantity 

r EE ffV^o • (9) 
Making use of the relation ||7J), we find that 

p=[{l + zy -3r]H^, (10) 

_^,2 r' - 3»o,„(l + z fl 

Px ^ , , (11) 



0,„ 



r(l + z) 
1 /' 



2r l + z 



5L 



"2"°'H ^W^j r ^^^^ 

where a prime represents a derivative with respect to z. 
Eliminating the 4>'^px term from Eqs. Hll|l and H12|l . we 
obtain 



4r(l + z)A 



5'm±^5'^-Sr{l + z)AB , (13) 
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where 



A. Case of p = f{X) - V{(l>) 



771-*- 

2r[r' -Snomil + z^I] 



X" 

o,„ - 



r 

2^ 



5L 



From Eq. H13|l we require the foUowing condition 
5[l > 8r(l + z)AB . 



(14) 



(15) 



(16) 



Once we know r and 6m in terms of z observation- 
aUy, Eq. (|13|l is integrated to give the functional form 
of I{z). It is worth mentioning that the function I{z) 
is determined without specifying the Lagrangian density 
p{(j),X). By using Eqs. H10() and (|11|1 . we obtain p and 
4>''^px as the functions of z. The energy density of the 



scalar field, p 
we find 



Px ~ P: 



is also known. From Eq. Q 



(In/)' 



(17) 



This case includes quintessence [/(^) = X\ and a 
phantom field [/(X) = -~X\. From Eq. |[TT|l we find 



/^3rio„(l + z)2/ 
r(l + 2) 



(20) 



If we specify the function f{X), the evolution of 4>'{z) 
and (/)(z) is known from r(z) and /(z). By Eq. (|17|l we 
find the coupling Q in terms of z and (p. Finally Eq. H10(l 
gives 



V = f+ [3r - (1 + z)r']i7i 



(21) 



The r.h.s. is determined as a function of z. Since z is 
expressed by the field 0, one can find the potential V{(j)) 
in terms of (j). In the case of quintessence without a cou- 
pling Q, this was carried out by a number of authors 
111 lia lia III lia- We have generalized this to a more 
general Lagrangian density p — f{X) — V{(f>) in the pres- 
ence of the coupfing Q. 



Hence the coupling Q is determined once / and (j)' are 
known. We need to specify the Lagrangian density 
p{(f>, X) to find the evolution of </)' and Q. 

The equation of state for dark energy, w = p/p, is given 

by 



B. 



Case of f{X)V{<j)) 



P 



(p^px - p 

(1 + z)r' - 3r 
3r - 3r!o,„(l + z fl 



(18) 
(19) 



A normal scalar field corresponds to w > —1, which 
translates into the condition px > by Eq. (|18|1 . From 
Eq. ((TT1) we find that this condition corresponds to r' > 
3riom(l + z)"^!, which can be also checked by Eq. ((T^ . 
Meanwhile a phantom field corresponds to px < and 
r' < 3riom(l + z)^/. As we already mentioned, the evo- 
lution of I{z) is determined if r and 5^ are known ob- 
servationally. Then the equation of state of dark en- 
ergy is obtained by Eq. H19|) without specifying the La- 
grangian density p{(j),X). In the next section we shall 
apply our formula to several different forms of scalar-field 
Lagrangian. 



III. APPLICATION TO SPECIFIC CASES 

Most of scalar-field dark energy models can be clas- 
sified into two classes: (A) p = f{X) — V{(j)) and (B) 
p = f{X)V{(j)). There are special cases in which cosmo- 
logical scaling solutions exist. This corresponds to the 
Lagrangian density (C) p = Xg{Xe^'*'), where A IS a con- 
stant and g is an arbitrary function. In what follows we 
shall consider these classes of models separately. 



This case includes k-essence and tachyon ^. The 
pressure density of the form 

p{X, ip) = K{ip)X + L{ip)X^ , X ^ ~[V^f/2 , (22) 

is transformed to the Lagrangian density p = f{X)V{(j)) 
with f{X) = -X + X^ and ^(0) = K^/L by field redef- 
initions: 



(23) 



We note that dilatonic ghost condensate model cor- 
responds to a choice K{ip) = — 1 and L{ip) = ce^"^ . In 
this case the potential of the scalar field has a dependence 
V{(j)) cx [Ulil by the field redefinitions 

For the Lagrangian density p = f{X)V{<f>) we obtain 
the following relation from Eqs. (|10|l and Hll() : 




,,2/x ^ /-3r!o„,(l + z)2/ 
/ r(l + 2)[(l + 2)r'-3r]i/2 ' 



V = 



[(1 + z)r' ~ 3r]gg 
/ 



(24) 
(25) 



Once we specify the form of f{X), one can determine the 
fimctions 4''{z) and (j){z) from Eq. (|24|l . Then we obtain 
the potential V{(f)) from Eq. ^E^. 



^ We note that a coupled dark energy scenario betw een k-essence 
field and dark matter was recently studied in Ref. |41| . 
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C. Case of p = Xg{Xe^''') 

When we construct realistic dark energy models, scal- 
ing solutions may play an important role to solve coin- 
cident problems. In this case the energy density of the 
scalar field is proportional to that of a barotropic fluid 
{p oc pm)- It was shown in Refs. [T3. l43| that the exis- 
tence of scaling solutions restricts the form of the scalar- 
field Lagrangian to be 



p^XgiXe 



A0N 



(26) 



We note that this property holds both for coupled and 
uncoupled models of dark energy. 

For the Lagrangian density (|26|l Eqs. (|10|1 and give 



(1 + zy - 3f7om(l + zfl 



g 2[(1 -I- z)r' - 3r] 

,2 _ 2[(1 + zy - 3r] 



r{l + z)^g 



(27) 
(28) 



where Y = Xe^'^. If we specify the functional form of 
g{Y), one can determine the function Y = Y{z) from 
Eq. (gH). Then we find </>'(z) and (/)(z) by Eq. The 
parameter A is known by the relation Y = {l/2)(p^e''"^. 

It was shown in Refs. |43 | that the quantity Y is 
constant along scaling solutions, in which case the l.h.s. 
of Eq. (|27|) is constant. In the absence of the coupling 
Q, i.e., / ~ 1, accelerated expansion is not realized for 
scaling solutions. This means that the r.h.s. of Eq. H27|) 
should not be constant for Q = when we use observa- 
tional data. If the coupling Q is present, there is a pos- 
sibility to find a situation in which the r.h.s. of Eq. 1271) 
does not vary. This corresponds to the case where accel- 
erated expansion is realized by scaling attractors. Thus 
one can directly check whether the present acceleration 
originates from scaling solutions with a non-zero coupling 
Q once we obtain accurate observational data of r{z) and 



IV. EXAMPLES OF RECONSTRUCTION 

In this section we show concrete examples of recon- 
struction for several dark energy models. We shall use 
the following parametrization for the Hubble parameter 



r{x) = nomx^ +Ao + Aix + A2X^ , 



(29) 



where x = 1 + z and = I — Ai — A2 — f^om- This 
corresponds to the following expansion for dark energy 



p = pqc {Aq + Aix + A2X^) , 



(30) 



where poc = 3i?Q. 

For a prior f^om = 0.3, the Gold dataset of SN obser- 
vations gives Ai = -4.16 ± 2.53 and A2 = 1.67 ± 1.03 



|47j . We note that the weak energy condition for dark 
energy, p> and w = p/ p > —1, corresponds to 



Ao + Aix + A2X^>0, Ai+2A2X>0. 



(31) 



If we use the best-fit values Ai = —4.16 and A2 = 1.67, 
for example, we find that the second condition in Eq. H31|l 
is violated around present {x ^ 1). This means that the 
field behaves as a phantom {w < —1). In the case of an 
ordinary scalar field such as quintessence, we need to put 
a prior Ai + 2A2X > 0. 

The reconstruction program of a quintessence-type 
scalar field has been already carried out in Ref . ^3 1 so 
we do not repeat it here. We shall study the cases of 
tachyon and generalized ghost condensate when the cou- 
pling Q is zero. We have not yet obtained good dataset 
about 6m{z), so the coupling Q is not well determined by 
current observations. 



A. Tachyon 

Tachyon corresponds to the case / 
Eqs. (|23|l and (j^SJ. Then we have 



-(l-2X)i/2 in 



r' - 3f7o™a;2/ 



V = 



H^rx[rx{r' - SflomX^I) + 3r - 
(3r - xr')H^ 



^1 - r^x^H^ct)'^ 



3T' (32) 



(33) 



We note that the reconstruction equations were derived 
in Ref. 0| for Q = 0. The equation of state for tachyon 
is w = 0^ — 1, which means that w > —1. Hence we 
should impose the prior given by Eq. (|31|l . 

In Fig. n we show one example of the reconstruction 
of the tachyon model for Q = (/ = 1). The field value 
is chosen to be (/) = at present (z = 0). The tachyon 
potential needs to be flat around < z < 1 to give rise to 
an accelerated expansion. The potential has a minimum 
with a non-zero energy density for the parametrisation 
given in this flgure {Ai = -3.80 and A2 = 1.95). This 
implies that the rolling massive scalar field model with 
potential V{(j)) = Vb exp(m^(/)^/2) can be a viable 
dark energy model. 



B. Generalised ghost condensate 

When the field satisfies the condition px > 0, we need 
to put the prior for consistency. Let us consider a 
situation in which crossing of the cosmological-constant 
boundary is possible. This can be realised for the follow- 
ing type of Lagrangian density: 



p= -X + h{(j>)X^ 



(34) 



where /i(0) is a function in terms of (f>. Dilatonic ghost 
condensate model corresponds to a choice h{(j)) = 



5 






FIG. 1: Reconstruction of tachyon model for Q = for the 
parametrisation 12911 with Ai = —3.80 and A2 — 1.95. These 
coefficients satisfy the weak energy condition (1311 . which 
means that the equation of state ranges in the region ui > — 1. 
We show V, w and z in terms of the function of We note 
that the potential is normalised by HQirn^i, where mpi is the 
Planck mass. 



FIG. 2: Reconstruction of generalized ghost condensate 
model for the parametrization 12911 with the best-fit parame- 
ters A\ = —4.16 and A2 = 1.67. We show h, w and z in terms 
of the function of (j>. This model allows a possibility to cross 
the cosmological-constant boundary {w = —1). 



ce^'^. From Eqs. and idJ we obtain 



h{<j>) 



12r - 3xr' - SriomX^I 



(35) 
(36) 



In Fig.|21wc plot h{(j)) in terms of the function when 
we use the best-fit values of ^1 and A2. The cross- 
ing of the cosmological-constant boundary corresponds 
to hX ^ 1/2, which occurs around the redshift z = 0.24 
for the best-fit parametrisation. The system can enter 
the phantom region [hX < 1/2) without discontinuous 
behaviour of h and X. 

However we have to caution that the perturbation of 
the field (p is plagued by a quantum instability whenever 
it behaves as a phantom [l2|. Even at the classical level 
the perturbation is unstable for 1/6 < hX < 1/2, since 
the speed of sound, = Px/{px + ^Xpxx), becomes 
negative. One may avoid this instability if the phantom 
behaviour is just transient. In fact transient phantom 
behavior was found in the case of dilatonic ghost conden- 
sate model (see, e.g.. Fig. 4 in Ref. ^3)- ^^^^ '^^^^ *he 
cosmological-constant boundary crossing occurs again in 
future, after which the perturbations become stable. 

We found that the function h{(j)) can be approximated 
by an exponential function e"*"^ near to the present, al- 
though some difference appears for z ^ 0.2. However 
the current observational dataset is not still sufficient to 



rule out the dilatonic ghost condensate model. We hope 
that future high-precision observations will determine the 
functional form of hi^cp) more accurately. 



C. Scaling solutions 

As we already mentioned, the existence of scaling so- 
lutions can be found by evaluating the r.h.s. of Eq. 127|) . 
When Q = we checked that the r.h.s. of Eq. l(77|l 
is not constant when we use the Gold dataset with the 
parametrisation (|29|l . which is consistent with the fact 
that scaling solutions do not lead to an accelerated ex- 
pansion for (5 = 0. 

In the presence of the coupling Q if the solution around 
< z < 1 corresponds to a scaling solution, the r.h.s. of 
Eq. (|77|l is constant. This gives the constraint for the 
evolution of /: 



liz) 



(37) 



where /(O) — 1. The r.h.s. of this equation is indepen- 
dent of the scalar-field Lagrangian. 

If both the evolution of r(z) and (Sm(z) are known, 
one can determine I(z) by Eq. H13|l . Then by comparing 
this with Eq. 1)3 7() , one can check the existence of scaling 
solutions. It was shown in Ref. [s^l that in the case of 
a phantom field the final attractor does not correspond 
to scaling solutions but to scalar-field dominant solutions 
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with fi^ = 1 . Hence we have to caution that Eq. (|37|l can 
be used for the region characterised by Eq. 



V. RELATIONSHIP WITH GENERALISED 
EINSTEIN THEORIES 

Having considered the couphng Q between dark en- 
ergy and dark matter, we would hke to relate this sce- 
nario with theories giving rise to such a coupling. Let us 
consider the following 4-dimensional Lagrangian density 
with a scalar field and a barotropic perfect fluid: 



1 



1 



(38) 



where F{ip), C{^p) and U{ip) are the functions in terms of 
if. This includes a wide variety of gravity models-such as 
Brans-Dicke theories, non-minimally coupled scalar fields 
and dilaton gravity^. 

Making a conformal transformation = F{ip)g^i,, 
the above action reduces to that of the Einstein frame 



^(V0)^ 



where 



^ / G(^)d^, G(^)^ 



(39) 



(40) 



We note that several quantities in Einstein frame are re- 
lated with those in string frame via a ~ 
VFdi, pm = Pm/F^, = pm/F^ and V = U/F^. In 
Einstein frame the background equations are given by 
Eqs. © and Q with p = X -V{(j)) and 




3H(p. 



2FG 



{Pm - 3Pm)4> ■ (41) 



In the case of non-relativistic dark matter (p„ 
corresponds to Eq. Q with a coupling 



F^ 
2F 



3 fF 
2 



F 



-1/2 



0) this 



(42) 



For example a nonminimally coupled scalar field with 
a couphng ^ corresponds to F{ip) = 1 — ^<p^ and ({ip) = 1. 
In this case we find 



(43) 



^ We note that the Lagrangian 1381 can be extended to the case in 
which hig her- order curvature corrections are taken into account. 
See Refs. l5ll for details. 



Then ^ is expressed in terms of Q: 



2(1-6Q2) 



-1 ± 



4(1 -6g2) 



(44) 



Once we know r(z) and 6miz), the coupling ^ is evaluated 
by using Eqs. (|17|l and (|44|l together with the relation 
(I40|l between and ip. Hence it is possible to determine 
the strength of the nonminimal coupling from observa- 
tions. We note that Q ~ ^ip for |^| <^ 1 and Q ~ ±1-^6 
in the large-coupling limit (|^| 3> 1). The large-coupling 
case is excluded from solar system experiments provided 
that the scalar field is universally coupled to all matter 
[23. [25I I . It is certainly of interest to place constraints on 
the strength of ^ by using our reconstruction formula to- 
gether with other experiments about the time-variation 
of a gravitational "constant" G. 

String theory also gives rise to the coupling Q after a 
conformal transformation from strin g fra me to Einstein 
frame. The tree-level dilaton gravity |53| corresponds to 
F{(p) ~ e~'^ and C(v) ^ "e""^! which gives a constant 
value of Q by Eq. (|42|) . It is typically assumed that non- 
perturbative effects would stabilize the dilaton with a 
potential so that it does not contradict with solar system 
experiments. An alternative possibility is the runaway 
dilaton scenario [s^l in which the dilaton is effectively 
decoupled from gravity in the limit (/? — > 00 with the de- 
pendence F{ip) - Bi + Cie-'fi and ((ip) ^ B2 + Cae^"^. 
One can also check the viability of this scenario by com- 
paring Eq. (|42|l with the coupling Q obtained from ob- 
servations. 

There is another interesting cosmological scenario in 
which neutrinos are coupled to dark energy [ssj . In this 
model the neutrino mass, m^, is a function of a scalar 
field, (j). In a situation where neutrinos are coUisionless, 
the neutrino energy density, p^, satisfies the equation of 
motion l55l 



p^ + 3H{p 



91n7 



— [pv - ipy 



(45) 



where p^, is the pressure density of neutrinos. When the 
neutrinos become non-relativistic, Eq. (|45|l shows that 
the coupling between neutrinos and dark energy is given 
by (5(0) = dhnriiy/dcj)- Hence if we know the coupling 
Q{4>) observationally, the evolution of the neutrino mass 
is found as a function of (f> (and z). 



VI. CONCLUSIONS 

In this paper we have provided a method to recon- 
struct scalar-field Lagrangian in an accelerating universe 
from observations. Our starting point is the general La- 
grangian which is the function of a scalar-field 4> and 
a kinematic term X. We have also taken into account the 
coupling Q{(f>) between dark energy and a non-relativistic 
perfect fluid in order to include the coupled quintessence 
scenario. 
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In the absence of the coupHng Q, one can reconstruct 
the structure of theory by parametrising the Hubble rate 
H in terms of the redshift z from the luminosity distance 
of supcrnovae observational data. We need to know ad- 
ditional information in order to determine the coupling 
Q from observations. We have made use of the equa- 
tion for matter density perturbations 5m on sub-Hubble 
scales for this purpose. Our reconstruction formula is 
given by Eqs. (|10(l . and H13() . from which the cou- 
pling Q is also determined with the use of Eq. H17(l . We 
note that the equation of of state l|19(l for dark energy 
is known from observables only without specifying any 
Lagrangian. 

In Sec. mil we have applied our reconstruction for- 
mula for several forms of Lagrangian density: (i) p = 
f{X) - (ii) p = f{X)V{<ty) and (iii) p = XgiXe^"^) 

where g is an arbitrary function. In the cases (i) and 
(ii) reconstruction equations can be decomposed into two 
contributions coming from a kinematic term and a poten- 
tial term. Hence the potentials of such theories can be 
obtained together with the coupling Q once the evolution 
of H{z) and 5m{z) is known. The case (iii) corresponds 
to the Lagrangian density for the existence of scaling so- 
lutions. One can check the existence of scaling solutions 
if we evaluate the r.h.s. of Eq. H27(l using observational 
data. 

In Sec. IIVI we have presented concrete examples of 
our reconstruction with the parametrization given by 
Eq. (|29|l . We studied two dark energy scenarios in the 
absence of the coupling Q- (a) tachyon and (b) gener- 
alised ghost condensate. In the case (a) the equation of 
state for the tachyon field is constrained to be w > — 1, 
which means that the prior needs to be imposed. In 
Fig. n we plotted one example for the reconstruction of 
the tachyon potential. The model (b) allows a possibility 
to cross the cosmological-constant boundary: w ~ —I. 



We carried out the reconstruction of this model by using 
the best- fit values coming from the Gold dataset. The 
result is illustrated in Fig. |21 which shows that the field 
behaves as a phantom for the redshift < z < 0.24. 

In Sec. we presented a Lagrangian in generalised 
Einstein theories which gives rise to the coupling Q by a 
conformal transformation to Einstein frame. For exam- 
ple, nonminimal coupling ^ is directly related with Q as 
given in Eq. (|44() . This allows a possibility to determine 
the strength of the nonminimal coupling from observa- 
tional data of H{z) and (5m(z). 

For the moment we have not yet obtained the accu- 
rate evolution of 5m{z) from observations of clustering. 
We only know the total amount of growth between the 
decoupling epoch and present. This is associated with 
the fact that all probes of clustering are plagued by a 
bias problem. However upcoming galaxy surveys such as 
KAOS, LSST and PANSTARS will pin down the matter 
power spectrum to exquisite accuracy, allowing the ulti- 
mate measurement of the power spectrum. By that time 
we should have an excellent understanding of bias and 
will be able to obtain the time-evolution of 5m- We hope 
that this will provide us an exciting possibility to reveal 
the origin of dark energy. 
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